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Abstract 

We extend earlier treatments of holographic superconductors by studying 
cases where operators of different dimension condense in both 2 + 1 and 3 + 1 
superconductors. We also compute a correlation length. We find surprising 
regularities in quantities such as ojg/Tc where ojg is the gap in the frequency 
dependent conductivity. In special cases, new bound states arise corresponding 
to vector normal modes of the dual near-extremal black holes. 



1 Introduction 



It has recently been shown that a simple Einstein-Maxwell theory coupled to a 
charged scalar provides a holographically dual description of a superconductor [Il|2]. 
The gravity theory depends on two parameters, the mass m and charge q of the 
complex scalar field. If q is large, the backreaction of the matter fields on the metric 
is small and can be neglected. This probe approximation was used in [2\ and will be 
adopted here as well. (For a discussion of the dependence on q, see [3].) Our main 
goal is to explore how the properties of the superconductor depend on the scalar mass 
m. This determines the dimension A of the operator that condenses. In addition, the 
previous studies have focussed on a four dimensional gravity theory which is dual to 
a 2 + 1 superconductor. This is appropriate for many high Tc materials in which the 
superconductivity is associated with two dimensional planes. Here we also study the 
five dimensional gravity theory dual to a 3 + 1 superconductor. (For other work on 
holographic superconductors, see |1]-[T2].) 

We consider a variety of masses between = and the Breitenlohner-Freedman 
bound < marking the boundary of stability for a scalar field in ArfS*!^ We 
find qualitatively the same behavior for all masses. Below a critical temperature, a 
charged black hole becomes unstable to developing scalar hair. In the dual theory, 
this corresponds to the formation of a charged condensate. By perturbing the black 
hole, one finds that the conductivity cr{uj) in the dual theory has a delta function 
at = 0, and a gap Ug. These are standard properties of s-wave superconductors. 
There are special features which arise when saturates the BF bound. We will see 
that at low temperature and general m^, cr(a;) has poles in the complex frequency 
plane whose locations depend on m^. When = Tn^p these poles move onto the 
real axis inside the gap. This produces delta function contributions to SHe[cr] at 
certain frequencies. This suggests that the charged quasiparticles are interacting 
with a strength that depends on w?. When w? = m\p, they form bound states]^ 
From the bulk standpoint, this corresponds to a (vector) quasinormal mode which 
becomes an actual normal mode at low temperature and m? = rn^Qp- 

In the cases we consider, the dimension A of the operator condensing changes by a 
factor of three for 2+1 superconductors and a factor of two for 3+1 superconductors. 
Some quantities change significantly such as the minimum energy A for charged 
excitations. In particular, the BCS relation uog = 2A which held in [2] is not true 

^We expect theories with > to also be dual to superconductors, but in this case there is a 
numerical instability at large radius making it difficult to find solutions with the correct asymptotic 
behavior. 

■^We thank M. Fisher for this interpretation. 
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in general. However we have noticed some striking regularities. In particular, in 
all cases where A > Xbf, ^g/Tc ~ 8. This holds to better than 10% for both the 
three and four dimensional superconductors, as well as various scalar masses. Since 
the corresponding BCS value is 3.5, this shows that the energy required to break 
apart the condensate is more than twice the weakly coupled value. The holographic 
superconductors are indeed strongly coupled. We also find regularities involving the 
superfluid density and a correlation length i^^.. 

In the next section, we describe the gravitational theories we will study and show 
that scalar hair forms at low temperatures in all cases. In section 3, we compute 
the conductivity in the dual theory by considering electromagnetic perturbations of 
the black holes. Section 4 contains a discussion of other physical quantities, such as 
the energy gap A, the superfluid density n^, and a correlation length The final 
section contains a short summary. In a brief appendix we calculate the conductivity 
in the normal phase analytically, finding agreement with the numerical results. 



2 Abelian Higgs Model 

We study an Abelian Higgs model in AdSd+i- 

S = j (f^^x^ (^^F^.F^'' - |V^ - iA^I' - m^l^'l^ . (1) 

(Note that we do not include a term.) If one multiplies this action by and 
sets \l/ = , one recovers the standard action for a scalar ^! with charge q. In the 
limit where q is large and \l/ is held fixed, the backreaction on the metric becomes 
negligible. Following [2], we will work in this probe limit. 

Our background metric will be the standard AdS-Schwarzschild black hole 

ds^ = -fdt^ + r^dr^ + r^dxidx\ f{r) = r\l - (2) 

where we have chosen units such that the AdS radius is unity. To see the formation 
of scalar hair, we consider static, translationally invariant solutions to ([T]). Consider 
the ansatz = {"^l = tp, A = (p dt where ip, are all functions of r only. This leads 
to the equations of motion 

r + —j ^' + y^V' - y = 0, (3a) 

0" + —0' -^0 = 0. (3b) 
r f 
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As pointed out in [T], the coupling of the scalar to the Maxwell field produces a 
negative effective mass for if). Since this term becomes more important at low tem- 
perature, we expect an instability toward forming scalar hair. This was indeed seen 
for = —2 and ci = 3 in [2]. We now consider = and = m%p = —9/4 
in d = 3 and = 0, —3, —4 in ci = 4. Note that the last value of is at the BF 
bound. Since ^ are coupled nonlinear equations, we do not expect to solve them 
analytically. However, it is straightforward to solve them numerically. 
Recall that the asymptotic behavior of a scalar field in AdSd+i is 

^ =— + — + •• • (4) 

where X± = \{d± ^/d"^ + Am?). For > — (i^/4 + 1, only the A4. mode is normal- 
izable, and so we interperet ip+ = {O), where O is the operator dual to the scalar 
field, and A is the dimension of the operator]^ is dual to a source for O. To 
have spontaneous symmetry breaking, we wish for the operator to condense without 
being sourced. For —d'^/4 < m? < —d'^/A + 1 both modes are normalizable and we 
have a choice of quantization. However, stability of AdS dictates that the only stable 
quantizations which respect AdS symmetries at large radius are ones where ■?/'+ or 
ip_ vanish [13]. When = —d^jA = m\p, A_|_ = A_ = Xbf = d/2 and there is a 
logarithmic branch. Such a branch will necessarily introduce an instability unless it 
is treated as the source [H]. We have not obtained the Abelian Higgs model from a 
string theory truncation, so the dual operator O is not known beyond the fact that 
it is charged under a global U{1) symmetry and has conformal dimension A at strong 
coupling. 

The asymptotic behavior of cf) is 

In terms of the dual theory, /i is the chemical potential and p is the charge density. 

In figure [l] we plot the condensate {O) for various masses and spacetime di- 
mension. The previously discussed case th? = —2 is included for comparison. It 
is convenient to label the curves by the dimension of the operator A rather than 
the mass m in order to distinguish the two possible quantizations when = — 2 
and d = 3. In all cases there is a critical temperature, Tc, above which there is no 
condensate. The critical temperature is proportional to p^/^ in d = 3 and p^^^ in 
d = 4, refiecting the different dimensions of the charge density in the two cases. The 
ratios T^/ p^/^ and Tc/ p^^^ decrease as A increases. This is expected since increasing 



^This normalization of O differs from that of [2 by a factor of y/2. 
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d=3 



d=4 




Figure 1: The condensate as a function of temperature. A is the dimension of the 
operator O, and d is the spacetime dimension of the superconductor. \bf = —3/2 
for d = 3 and Xbf = —2 for d = 4. The condensate tends to increase with A. 



A corresponds to increasing the mass of the bulk scalar, and hence making it harder 
for the scalar hair to form. As expected from mean field theory, (O) ~ (Tc — T)^/^ 
near the critical point. At low temperature, for A > Xbf, the condensate quickly 
saturates a fixed value which increases with A. For A = Xbf, the vev approaches a 
fixed value roughly linearly, and for A < Xbf it appears to diverge. 



3 Conductivity 

To observe that our boundary theory is superconducting, we need to calculate the 
conductivity a. This is related to the retarded current-current two-point function for 
our global U{1) symmetry, <j{uj) = j^G^iu, k = 0). To do so we must calculate an 
electromagnetic perturbation on top of the hairy black hole. The linearized equation 
of motion for 5 A = Ax{r)e~'''^^^^^'^ dx is 

ff d-3\ ^, (U? e 2tP^\ , 

K,[^.—)K,[j,-^~^)a^O (6) 

This mode is not coupled to other hnearized perturbations and can be studied sepa- 
rately. We require oc f-^^/'^'^o near r = Tq corresponding to ingoing wave boundary 
conditions at the horizon. The AdS/CFT dictionary tells us how to calculate the 
retarded current Greens function from our gauge field perturbation [TJ]. The 
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result is 

= - lim fr^-'A,A', (7) 

r^oo 

where is purely infalling at the horizon, and normalized to Ax{r = oo) = 1. 

Let us work this out for the two cases of interest to us, c? = 3 and 4. As can be 
verified by the large r expansion of (|6|), in (i = 3, the gauge field falls off as 



4(1) 

A. = A(°) + — + ••• (8) 
r 

giving 

= ''(-) = ^''>-* = «' = ^i'- 

The d = 4 case is a little more nuanced. There is a logarithmic term in the falloff 
at nonzero frequency and momentum, 

^, = ^<°» + ^+-^'°'<-'-^')'ggj^... (10) 

giving 

A (2) 1 

G^ = 2j^ + {u'-e){logAr--). (11) 

This logarithmic divergence can be removed with a boundary counterterm in the 
gravity action [16], but it necessarily breaks conformal invariance. In other words, 
we must specify a renormalization scale when regulating the action. After adding a 
counterterm to cancel the logAr term we have the Greens function 

It is important to note that when we are calculating physical quantities such as 
Un, rig, A, Ug, these will all be invariant under a shift in A, as all this can change is 
the quadratic term in 9^e[G'^]. It will however effect the definition of our correlation 
length, as will be explained later. 

In Figures [2] and [3] we plot the conductivity o-{tu) for the various cases studied. It 
is clear that in all cases we find a gap in the conductivity. We parameterize this by the 
gap frequency Ug. At zero temperature we expect there to be a sharp gap, indicating 
that there are no states below said energy scale. In most cases the conductivity rises 
quickly near Ug, the exception being the case when the bulk scalar is massless. In 
this case, while there is exponential suppression in cr at -C ujg, the conductivity 
rises more gradually around the gap frequency even at low temperature. 
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Figure 2: Conductivity for 2 + 1 dimensional superconductors. Each plot is at low- 
temperatures, about T/Tc ~ 0.1. The solid line is the real part, dashed is imaginary. 
The pole at = is clearly visible in Jm[(j]. For A > 3/2, we find the shape of the 
curve near the edge of the gap largely dictated by a second pole in the lower half 
complex io plane. As T — > and A — > 3/2, this pole hits the real axis. 
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Figure 3: Conductivity for 3 + 1 dimensional superconductors. Each plot is at low- 
temperatures, about T/Tc ~ 0.1. The solid line is the real part, dashed is imaginary. 
The pole at cij = is clearly visible in Dfm[(T]. For A = 2, we find the shape of the 
curve largely dictated by multiple poles. As T ^ and A — > 2, these poles approach 
and hit the real axis. 
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At nonzero temperature, there is some ambiguity in defining Ug when looking at 
the thermally smoothed conductivity plots. However, it turns out that for A > Xbf, 
Jm[cr] always has a minimum right in the region where 9^e[cr] is turning on. We 
will define Ug as the frequency minimizing Jm[o"]. For A < Xbf, 3ftn[(T] decreases 
monotonically and we define Ug by the minimum of \a\. 

We find that ujg remains quite uniform for many of the cases we have studied. 
For all cases with A > Xbf, we find 




with deviations of less than 8%. Even changing the spacetime dimension of the su- 
perconductor from d = 3 to d = A does not affect this ratio. Since the corresponding 
BCS value is 3.5, this shows that the energy to break apart the condensate is more 
than twice the weakly coupled value. 

The most interesting case is A = Xbf- Here we find that additional poles in 
Jm[(j] and delta functions in ^Rc[a] arise at low temperature inside the gap. One can 
interpret such poles as follows. In standard BCS theory, the quasiparticles excited 
above the superconducting ground state are not interacting. Since holographic su- 
perconductors are strongly coupled, it is perhaps not surprising to find interactions 
between the quasiparticles. For A > Xbf, the Greens function has poles in the 
complex frequency plane which produce the bump in 9^e[cr] just above the gap. This 
looks like a resonance in quasiparticle scattering. As A ^ A^f and T ^ 0, the pole 
moves onto the real axis, indicating that the interactions become stronger and form 
a bound state. For d = 3, there is one such pole, and the difference Ug —Upoie reflects 
the binding energy of the bound state. For d = 4, more poles appear as one lowers 
the temperature. Since our numerics can only explore down to T/Tc ~ 0.1 it is not 
clear what the total number is, and what the ultimate gap Ug will be. It would be 
very interesting to study these cases with backreaction to see how much of this pole 
structure remains. 

On the gravity side these poles are also surprising, as they imply that our near 
extremal black hole has normal modes. More precisely, at the frequency Upoie there 
are solutions for a vector potential which are ingoing at the horizon and vanish at 
infinity. Usually, black holes have quasinormal modes with complex frequencies that 
describe the decay of perturbations outside the horizon. In AdS, these correspond 
to poles in the retarded Greens function flbl [T7] . Remarkably, when = m'^p and 
T ^ 0, some quasinormal frequencies become real. There are vector perturbations 
which do not decay. We do not have an intuitive explanation for this. It would be 
interesting to understand why this is happening and check if there are similar normal 
modes associated with the scalar field ip as well. 
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Figure 4: The absolute value of the retarded Greens function is shown as a function 
of temperature and (real) frequency for A = A^f at zero spacial momentum. Lighter 
shades denote larger values. One can clearly see the complex frequency poles as they 
move onto the real axis. 



Figure |4] shows \G^\ as a function of temperature and (real) frequency for the 
A = Afii? cases. Lighter shades denote larger values. One can clearly see the poles 
as they move onto the real uj axis. One can infer the conductivity at different 
temperatures by looking at horizontal lines in this diagram. 

4 Other properties of the superconductor 

In BCS theory, one often studies the order parameter A(T), the mass of charged 
quasiparticle excitations about the fermi surface as a function of temperature. How- 
ever, the gravity dual does not have direct access to the elementary charged particles 
in the dual superconductor, so we cannot calculate this in general. We can, however, 
extract the mass of the low lying excitations at low temperature. We look at the 
normal contribution to the DC conductivity, n„ = lim^^o which will be ex- 

ponentially suppressed when these modes are gapped, n„ ~ ^-^i^)/'^ , This gives us 
the capability to numerically calculate A(0). Note that for (i = 4 the conductivity is 
dimensionful, and we read off A from a /T^. In the BCS theory, uOg = 2A, due to the 
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Figure 5: Superfluid density in d = 3 and d = 4. 



fact that the electron-photon couphng imphes that a photon will excite quasiparticles 
in pairs. We find that this BCS result does not hold in general, even though it did 
hold for the cases studied in 0. The values of A are given in Tables 1 and 2 at the 
end of this paper. 

Another order parameter is the superfluid density, Ug. It is the coefficient of 
the pole in Jm[o"] at = 0. By the Kramers-Kronig relations, if Jm[cr] = ris/uj, 
9^e[(T] = 7ms6{uj). It is this delta function that tells us our boundary theory exhibits 
DC superconductivity. The BCS theory gives relations between Us, (O), and A, but 
for holographic superconductors, no analogous formulas are known. 

It is clear from the deflnition of a that = lHe[G^(ti; = k = 0)]. In terms of 
long-wavelength response, we can extract two correlation lengths, and ^k- We 
simply study the small u, k behavior of G^, 

9^e[G^] = n,(l-e>2 + e^fc2 + ---) (14) 

Because Lorentz invariance is broken at flnite temperature, we do not expect C,lo = 
^k- However, at low temperatures the spacetime is in some sense "nearly" AdS in 
Poincare coordinates, and so at low temperatures we expect the two values to be 
approximately equal. In the d = 3 case we indeed flnd that they are nearly equal not 
only at low temperature, but all the way up to T/Tc ~ 0.7. For c? = 4, and are 
not separately well deflned since there is a renormalization ambiguity as mentioned 
earlier. We can always shift G^ + G{uj'^ — k"^). Therefore in that case the only 

renormalization-independent quantity is the difference, 5^"^ = — ^5- Studying G^ 
at small k can be understood as linear response to a long wavelength magnetic fleld, 
it is unclear how to interperet ^i^. 
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The results for the superfluid density and correlation length are given in Fig- 
ures 5 and 6. We find a surprising relation between Ug and the gap frequency ojg. 
For all A > \bf we find 

--^ (rf = 3), {d = A) (15) 

The agreement is better than 10%. For = 3, A = 1 we find that Us/ujg = 1. 
This agrees with the analytic zero temperature approximate solution in [3j. With 
ijjir) ^ (Oi)/r, one finds G^(T = 0) = v/2(Ci)2 - + P. Note that this also gives 
the prediction = = 73^- indeed find from numerics that ik{Oi) = 0.51. 
We also find that, for all four values of A studied in c? = 3, 

efe(0)(OA)'/^~.18 + .37A . (16) 

It would be very interesting to see how much of this universal behavior remains when 
backreaction is included. As mentioned above, we cannot calculate ^ki^x)^^^ in 
d = 4, and due to restoration of Lorentz invariance 6C, vanishes at low temperatures. 

We find universal behavior near the critical point. As was pointed out earlier, the 
condensate vanishes as (Tc— T)^/^, as expected from mean field theory. The superfluid 
density always vanishes linearly, as expected from Landau-Ginzburg theory, and the 
coefficient is roughly the same (within 15%) for each dimension: 

^20(Te-T) (rf = 3), ^ 100T,(T, - T) (rf = 4) . (17) 

Similarly, the correlation length diverges as (T^ — T)~^/^ as is expected from Landau- 
Ginzburg with nearly exactly the same coefficient in each dimension: 

a = ^(l-T/T,)-i/2 id = 3), a=°|^(l-T/T,)-i/2 (d = 4). (18) 

A correlation length defined from perturbations in the ip field was found in the 
A = 2 case to have the same divergence in [7], and one would expect it to generalize 
to arbitrary A. The difference is that the correlation length we calculate is simi- 
lar to that defined in BCS theory, i.e., the small spatial momentum correction to 
the current-current correlator, whereas Maeda and Okamura calculate the Landau- 
Ginzburg correlation length, obtained from fluctuations of the charged scalar field. 



5 Summary 

We have studied holographic superconductors in both 2 + 1 and 3 + 1 dimensions, 
with operators of different conformal weight condensing. One advantage of going to 
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Figure 6: Correlation length in d = 3 and d = 4. 



3 + 1 dimensions is that there are various subtleties in 2 + 1 dimensions associated 
with the Mermin- Wagner theorem stating that one cannot break continuous sym- 
metries and the existence of Kosterlitz-Thouless transitions. The fact that we see 
qualitatively the same behavior in both cases shows that these subtleties do not play 
any role here. This is probably because we are using classical gravity to describe the 
superconductor, which is valid in a large limit. 

A summary of our numerical results are enumerated in Table 1 and Table 2 below. 
There is no entry for Ug in the X = \bf cases because, as we have seen, additional 
poles arise as T — >■ and Ug increases. At the smallest T/Tc that we can explore, the 
conductivity has not reached a limiting curve. The value of A seems anomalously 
low for d = 3, X = 3. This is a case where the conductivity rises slowly near Ug even 
at low T and it is possible that this has contaminated our measurement of A. 

We have found a number of surprising regularities involving the gap frequency 
ujg, the critical temperature Tc, the superfluid density Ug and the correlation length 



^k, see eq. (13) and (15) - (18). Some of these are clearly visible in the tables, such as 
the last line, and the row for ujg. Others involve comparing two quantities. It would 
be interesting to investigate this further and determine if some of these regularities 
lead to universal predictions for holographic superconductors. 
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Table 1: Various 2 + 1 superconductors. There is no entry for when A = 3/2 due 
to comphcations arising from complex poles moving onto the real axis (see text). 
Note t = T/Tc- 
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Table 2: Various 3+1 superconductors. There is no entry for ujg when A = 2 due to 
complications arising from complex poles moving onto the real axis (see text). Note 
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A Normal phase solution 

In this appendix, we calculate the conductivity in the normal phase, when the con- 
densate vanishes, for both d = 3 and d = 4. It turns out that this can be done 
analytically. At zero spatial momentum with no scalar hair, one can solve the wave 
equation for exactly. The infalling solutions are 



AJr) 



(( 



exp 



6ro 



log 



[r - To) 



+ rro + Tq 
1 



arctan 



4 4 ro 4 4 ro 



2r + ro 
V3ro 

ioj 



{d = 3) (19) 



Expanding at large r yields the normal phase conductivities 

d=3 



a. 



47rT 



+ % 



1 



:i+zv 

47rT 



, 2V27rT ■ 
+ log — — h 7 



V2](rf = 4) 
(20) 

(21) 
(22) 



-ttT, 



where here '?/'(z) is the digamma function '^{z) = T'{z)/T{z) 
For < T, 



a. 



T vr + z lo£ 



\ Tc ) f 



24T2 



(23) 



and for ^ T, 



d=4: 



VT . UJ 

2+^(^°^27; 



7) 



+ 



(24) 



It is clear from (19) that this particular perturbation has no homogenous quasi- 
normal modes in AdS^, though they exist at finite spatial momentum, as observed 
in [18]. This is due to the fact that the wave equation for a massless p-form field 
in d + 1 dimensional planar AdS-Schwarzschild reduces to a one- dimensional wave 
equation which has a fiat potential when d = 1 + 2p and k = 0. In particular, (19) 
can be written A^(r) = e"*"''* where r,, = J dr/f. 
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The AdS^ case has quasinormal modes at 

^QNM = 27rT(l - i)n, n = 1, 2, ■ • • (25) 

As was noted in [T7]. 

Our numerical results agree with these analytic formula for T > Tc, as well as 
T < Tc and u ^ Ug. This is expected, as at large u the frequency term dominates 
over the scalar hair contribution in the wave equation (Ig]). 
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